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I review the extraction of kinetic and chemical freeze-out parameters from experimental 
data, with particular emphasis on the underlying assumptions and the validity of the 
conclusions. 

1. INTRODUCTION 

Collisions of elementary particles, hadrons, and nuclei at ultrarelativistic energies pro- 
duce a multitude of particles ( "multiparticle production"). Final-state interactions be- 
tween the produced particles determine the dynamical evolution of the system. In e + e~ 
and hadron-hadron collisions only few particles are produced, and it is unlikely that many 
final-state interactions occur. The particles decouple ( "freeze out" ) from the system soon 
after production. 

On the other hand, in AA collisions the density of produced particles is sufficiently 
large over an extended region in space-time, such that the mean free path of produced 
particles becomes small and many final-state interactions occur. These interactions drive 
the system towards local thermodynamic, i.e., thermal, mechanical, and chemical equilib- 
rium. In local thermodynamic equilibrium, the evolution of the system is governed by the 
equations of ideal fluid dynamics Jl]. If the system is only thermally and mechanically, 
i.e., kinetically, but not chemically equilibrated, these equations have to be supplemented 
by rate equations which determine the chemical composition of the system 0. In both 
cases, pressure gradients between dense, equilibrated matter and the vacuum drive collec- 
tive expansion, which cools and dilutes the system. Freeze-out of particles occurs when 
microscopic interaction rates become smaller than the macroscopic expansion rate of the 
system. 

By definition, after freeze-out the momenta of the produced particles do not change. 
The experimentally measured spectra of hadronic particles thus reflect the state of the 
system at freeze-out. The question is whether the spectra can also tell us about the state 
of the system prior to freeze-out? For instance, can they tell us whether the system was 
in thermodynamic, or at least kinetic equilibrium? Do they provide information as to 
whether a quark-gluon plasma (QGP), i.e., an equilibrated state of quarks and gluons, 
was created at some stage during the evolution of the system? 



*Permanent address since February 2001: Institut fur Theoretische Physik, Johann Wolfgang Goethe - 
Universitat Frankfurt, Robert-Mayer-Str. 8-10, D-60054 Frankfurt/Main, Germany 



2 



These questions are addressed in the following. In section |, I discuss whether the 
exponential nature of single-inclusive, invariant particle spectra is sufficient to conclude 
that the system was in thermal equilibrium. In section || it is argued that in AA collisions 
collective expansion is likely to occur prior to kinetic freeze-out. In section |], I discuss 
how final-state particle ratios yield information about the thermodynamic conditions at 
chemical freeze-out. Conclusions are given in section ^[ 

2. EXPONENTIAL PARTICLE SPECTRA 

In multiparticle production processes, single-inclusive, invariant particle spectra are 
typically exponential in the transverse momentum, px- The reason for the exponential 
behavior is the phase space of the final multiparticle state ||. Consider for instance a 
process where two particles (or nuclei) collide and produce N particles in the final state. 
Label the 4-momenta of the two incoming particles as p'^" and p' 2 M , and the 4-momenta 
of the N outgoing particles as p%, p%, . . . , jf N . All 4-momenta are on-shell, pf = (Ei, pj, 
E,i = Jpf + mf. The total 4-momentum is conserved, P^ 1 = p[^ + p' 2 ^ = Y^LiVi- Up to 



constants depending on p'^ and p' 2 ^, the total cross section for this process is || 



N / N 



°~ /n^ (4) ( pM -£^] \m(p, p1 ,..., pn )\ 2 , a) 



where the integration is over the 3 N dimensional momentum space of particles in the 
final state, the delta function represents energy conservation, and \M.\ 2 is the modulus 
squared of the matrix element for the 2 — * N process. The single-inclusive, invariant 
cross section for the production of a certain particle species - without loss of generality 
assumed to have 4-momentum p^ - is then (again, up to constants) 

An r N H 3 tv / N \ 

^^-/n^^f^-rf-Erf) \m(p, Pi ,..., pn )\ 2 . (2) 

^- % — 2 \ % — 2 / 

Let us introduce the "average" matrix element, 

2\ 



P2, — ,PN 



\M(P, Pl ,..., PN )\ 2 ) i 

r N d 3 n- / N \ 

= [HP-Pi)]- 1 li^^lP'-Pl-EPi) \M{P,Px, 



,Pn)\ 2 , (3) 



where the average is over the 3 (N — 1) dimensional momentum space of the unobserved 
N — 1 particles in the single- inclusive production of the particle with momentum pi, and 



r N H 3 n- / N \ 

hp -pi) = /n ^ ^ (4) [ p "-p" - ix) 



(4) 



is the corresponding Lorentz-invariant momentum-space volume. The momentum-space 
volume has dimension MeV 2 ^" 3 -* . Then, the single- inclusive cross section can be written 

as 

Ey ^- ~ $(P - Pl ) (\M(P, Pl ,...,p N )\ 2 ) . (5) 

d d px ^ / p 2 ,...,PN 
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Dynamical information about the 2 — > N process is contained in the second term only. 
The first factor, the momentum-space volume, is in some sense trivial. 

To compute the momentum-space volume, for the sake of simplicity (and because it 
allows us to obtain a purely analytic result) let us assume that all N particles in the 
final state are massless. In this case, no other scale with the dimension of energy enters 
$(P — pi), and for dimensional reasons, 

r n N— 3 

$(P- Pl )~ [(P-Pi) 2 } ■ (6) 

Because $(P — p\) is Lorentz-invariant, one may evaluate the right-hand side in the CM. 
frame of the incoming particles, where P M = (E, 0), with the result 

*(P - Pi) ~ P 2(7V - 3) ( 1 - ) ■ (7) 

V F1J \ E/2N N N-3J K J 

In the limit N 3> 1, the term in parentheses is a representation for the exponential 
function, lim n ^ 00 (l + x/n) n = e z , such that 

$(P - ft) ~ exp , N » 1 . (8) 

Now insert this into Eq. ([|), and divide by the total cross section to obtain the invariant 
momentum spectrum. Then, for a given rapidity yi, for instance y\ — 0, one obtains 



diV 



dyi d 2 p T; i 



yi=o 



Provided that the average matrix element squared (|A^| 2 ) does not contain a strong (ex- 
ponential) dependence on p^], the invariant momentum spectrum decreases exponentially 
with px,i, with an inverse slope parameter T s i ope = E/2N, which proves our original as- 
sertion. The exponential behavior of the transverse momentum spectra is due to phase 
space 0. No assumption about thermal equilibration is necessary to obtain this result. 

Nevertheless, single-inclusive particle spectra are also exponential in thermal equilib- 
rium, i.e., for a system at temperature T. In this case, the invariant transverse momentum 
spectrum is given by the Cooper-Frye formula 

E 1 gL=ftt.p l f(a^,\ l ), do) 



d 3 p! 7s V T 

where S is the 3-dimensional space-time hypersurface on which the transverse momen- 
tum spectrum is computed, dS M is the normal vector on E, and f(x, A) is the thermal 
distribution function. For Boltzmann particles, f(x,X) ~ A exp(— x); A = exp(/z/T) is 
the fugacity of the particle, and fi the chemical potential. The 4-vector in Eq. ( [3~Q|) is 
the 4- velocity of the system, i.e., the average 4- velocity of particle flow. In the rest frame 
of the system, m m = (1, 0). 

For the sake of simplicity, compute the spectrum at constant time, where dS M = 
(d 3 x, 0), such that 

Pi-^-^AxPx exp(-^-) . (11) 



d 3 Pl x x V T 
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For ultrarelativistic particles, the energy per particle is related to the temperature via 
E/N = 3T, such that at midrapidity 



This is rather similar to Eq. (§), except that the (inverse) slope of the p^-spectrum is equal 
to the true thermodynamic temperature T = E/3N, while previously, T s i ope = E/2N is a 
factor 3/2 larger. 

The origin of this discrepancy is that invariant momentum space is not identical to 
thermodynamic momentum space. For a single particle, the former is d 3 p/E, while the 
latter is d 3 p. The missing factor of 1/E is responsible for the difference in the exponential 
slope. 

From an experimental point of view, without a fully exclusive measurement one cannot 
precisely tell the number of particles in the final state, such that the energy per particle 
E/N is unknown. The important point is that then there is no possibility to distinguish 
between the two cases Eqs. and (|T^): usually, the strategy is to perform an exponential 
fit to transverse momentum spectra with a slope parameter T s i ope , but in this way one 
cannot test whether the relationship between E/N and T s i ope is the same as in thermal 
equilibrium. Although the spectra are exponential, like for a thermal system, the system 
need not be in thermal equilibrium. Although the slope parameter is a quantity with 
the same dimension as temperature, temperature is not defined, if the system is not in 
thermal equilibrium. One may even go one step further: although single-inclusive spectra 
for different particle species may have the same slope parameter, this does not mean 
that the system is in thermal equilibrium. It is simply due to the fact that the slope is 
proportional to the energy per particle E/N, which is a constant for all iV particles in 
the final state of a 2 — > iV process at a given CM. energy E. 

3. KINETIC FREEZE-OUT 

Slope parameters for the most abundant particle species show different behavior as a 
function of particle mass in AA as compared to pp collisions [0. In AA collisions, they 
increase linearly with the particle mass, 

T s iope, - a + b m , (13) 
while in pp collisions, they are independent of the particle mass, 

Tsiope, pp — O . (1^) 

For CERN-SPS energies, a ~ 140 MeV for both AA and pp collisions, while b depends on 
the system size; it increases with the mass number A of the nuclei. 

Although there might be other explanations for this behavior, the most natural in- 
terpretation is in terms of collective motion. The parameter a is the slope parameter 
determined by multiparticle momentum space. Since the energy per particle E/N is an 
intensive quantity, in both cases a ~ E/N is independent of the system size. Once the 
particles are created, however, the environment matters, i.e., whether there are many 




(12) 
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final-state interactions, like in AA collisions, or whether the particles more or less freely 
stream towards the detector, like in e + e~ and pp collisions. As explained in the introduc- 
tion, in the first case local kinetic equilibration of the system leads to pressure gradients, 
which in turn generate collective motion. The constant b parametrizes this collective mo- 
tion; it is proportional to the (average) collective flow velocity of the expanding system 
J7J. The larger the system, the more final-state interactions occur, and the larger is the 
collective flow velocity. Consequently, b increases with the system size. 

A word of caution is in order. The fact that there is collective motion in AA collisions 
does not necessarily mean that local kinetic equilibrium is established and, consequently, 
ideal fluid dynamics applies to determine the evolution of the system. Pressure gradients 
which drive collective motion occur also in systems away from equilibrium. Nevertheless, 
the fact that there is collective motion indicates that there are final-state interactions, 
which will eventually drive the system towards kinetic equilibrium, unless the macroscopic 
expansion rate considerably exceeds the microscopic scattering rate. Therefore, the ideal- 
fluid approximation for the dynamical evolution of the system might not be too bad, 
especially for collisions of large nuclei at ultrarelativistic energies. 

Freeze-out of particles occurs when microscopic interaction rates become small as com- 
pared to the macroscopic collective expansion rate. In fluid dynamical models, one com- 
monly assumes kinetic freeze-out to happen instantaneously along space-time hypersur- 
faces of constant density or temperature. This is certainly an idealization: microscopically, 
a particle has a certain probability to decouple anywhere and anytime during the evolution 
of the system ||. 

The invariant momentum spectra of frozen-out particles are computed according to 
Eq. flTOP, with the freeze-out hypersurface X = Sf . There are a number of conceptual 
problems with this formula, a discussion of which is beyond the scope of this talk, see 
Refs. || for more details. In practice, however, the Cooper-Frye formula fllPp is sufficient 
to compute the particle spectra to reasonable accuracy. 

Let us assume that kinetic freeze-out happens across a surface of constant temperature 
T f o . Applying the mean- value theorem to Eq. (|T0|), the spectrum of frozen-out particles 
at midrapidity y = is 

_ (7) m r -p,. ( v T ) j (15) 

where (vy) and (7) are suitably defined average values for the transverse fluid 3-velocity 
and the Lorentz gamma factor along the freeze-out hypersurface. At y = 0, it is reasonable 
to neglect longitudinal collective motion, such that 7 ~ (1 — vt 2 ) -1 ^ 2 is determined by 
vy. Assuming azimuthal symmetry, the spectrum ( |I5|) then depends only on two param- 
eters, the kinetic freeze-out temperature Tf. Q and the modulus of the average transverse 
collective flow velocity (i>t)- 

Table p] shows values for Tf. Q . and (vt) extracted by several authors from experimental 
spectra for PbPb collisions at CERN-SPS energies, = 158 AGeV. From Eq. ( |T5| ) it 
is obvious that Tf and (vt) are correlated: fits of similar quality can be obtained by 
trading off a lower value for Tf against a higher value for (v?) and vice versa. This 
can also be seen in the values quoted in Table |l|. This ambiguity can be removed by 
two-particle correlations, where Tf. Q . and (vt) are correlated in the opposite way fill . 



dN 



dy d 2 p 7 



cxp 
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Table 1 

Kinetic freeze-out temperatures and average transverse flow velocities for PbPb collisions 
at E Lah = 158 AGeV. 

Ref. Tf. . [MeV] (v T ) 



140 0.55 

[[TJ 110 - 120 0.60 

[0 120 - 140 

@ 130 0.50 

[|TJ 120 0.43 

[|TJ 100 0.55 

140 0.55 



4. CHEMICAL FREEZE-OUT 

In principle, there is another, independent way to determine thermodynamic quantities 
at freeze-out. First note that, in order to compute the transverse momentum spectrum 
fllBD of particle species i, one not only needs to know the temperature and the transverse 
flow velocity along the freeze-out hypersurface, which determine the shape of the spectrum 
as a function of pr, but also the fugacity A^f.o. = exp(/Xj i f. ./Tf .), which determines the 
absolute normalization of the spectrum, cf. Eq. ([[(]). In general, /z^f.o. varies along the 
freeze-out hypersurface, but in order to proceed, let us make the assumption (the first 
of three) that it is constant, like the temperature Tf. .. In the following, we shall drop 
the subscript "f.o.", but remember that all thermodynamic quantities, as well as the flow 
4-velocity w M , are taken on the freeze-out hypersurface. 

The second assumption we shall make is that all particle species freeze out across the 
same freeze-out hypersurface. This is certainly an idealization, as the mean free paths of 
different particle species are different |§. Under these two assumptions, the ratio of the 
total particle numbers of particle species i and j is 

N t _ J(dW^)/sdS-^/fa^/T,A t ) _ / E dS-M 

/(d'p,.. /•.' ) j, <1>; • Pj f(pj ■ u/T, Xj) j, .IX -A} • } 

where 

d 3 Pi uj. fPi-u 



AO 



Ei 



Pf/(^,A,) (17) 



is the 4-current of particle species i. In kinetic equilibrium the 4-current assumes the 
simple form Af- 1 = riiU^ 1 , where rij = rii(T, [ij) is the density of particles of species i in 
the local rest frame. The flow 4-velocity u^ 1 is common to all particle species because the 
system is assumed to be kinetically equilibrated immediately prior to freeze-out. Equation 



Q16D becomes |T 



Nj _ ni(T,fii) 
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because the factor / s dS • u cancels between numerator and denominator. The result 
fll8| ) is remarkable in the sense that, under the present assumptions, the ratio of total 
particle yields is independent of the detailed dynamical evolution of the system prior to 
freeze-out. Moreover, this ratio is the same as for a system in global kinetic equilibrium 

at temperature T with chemical potentials /ii, ^2, 

Let us now make the third assumption, namely that the system is not only in kinetic, 
i.e., thermal and mechanical equilibrium, but also in chemical equilibrium. This assump- 
tion is fulfilled if inelastic (and not only elastic) collision rates are sufficiently large. In 
this case, the chemical potential of particle species i is given by 

/ii = h fi B + Si /j s + e-i He + ■ ■ • , (19) 

where /zg, /is, and /i e are baryon, strangeness, and electric charge chemical potentials, 
respectively, and b iy Si, are baryon, strangeness, and electric charges of particle species 
i. The ellipsis in Eq. (|T^) stands for possible other conserved charges with associated 
chemical potentials. 

Global baryon, strangeness, and charge conservation impose additional conditions which 
allow to eliminate all charge chemical potentials except for one, e.g. fis- Then, all 
particle ratios are a function of two parameters only, T and [is- Taking into account the 
known (vacuum) values for the hadronic masses and the (vacuum) partial decay widths 
Ti_>jk..., one can then extract T and /i# from (47r extrapolated) data via a x 2 analysis. 

Different authors use slightly different strategies to perform such x 2 fits. Some introduce 



a hard-core repulsion between hadrons [|T^,|T^], some relax the assumption of chemical 
equilibration of strangeness [T2^,^T|], some conserve strangeness exactly in the canonical 
ensemble ||22|| . In essence, all methods introduce at least one additional parameter, which 



of course improves the quality of the fit, but does not fundamentally change the underlying 
assumptions. All fits are roughly of the same quality, which makes it hard to draw definite 
conclusions regarding the necessity of the individual approach. 

Results of x 2 analyses for e + e~, pp, as well as AA collisions with equal- mass nuclei at 
various collision energies are compiled in Fig. |l] (see also [17];|23| for a similar compilation). 
Two features are remarkable: first, as noted by Cleymans and Redlich |fT7j| , all (T, //g) 
combinations fall on a line of constant energy per particle, (E) / (N) ~ 1 GeV. The value 
1 GeV can be intuitively understood as setting the energy scale below which inelastic 
collisions cease and therefore chemical equilibration becomes impossible. 

The second feature is that, at CERN-SPS energies and above, the freeze-out tempera- 
ture is T ~ 160 MeV, which is somewhat higher than the kinetic freeze-out temperature 
discussed in the last section. This can be understood rather naturally noting that the 
extraction of a freeze-out temperature from particle ratios assumed chemical equilib- 
rium, while the extraction from single-inclusive spectra only assumed kinetic equilibrium. 
Thus, particle ratios determine the temperature at chemical freeze-out, and not at kinetic 
freeze-out. Since chemical equilibrium requires frequent inelastic collisions, while kinetic 
equilibrium only requires frequent elastic collisions, chemical freeze-out occurs earlier, 
i.e., at larger temperatures, than kinetic freeze-out. In order to distinguish quantities at 
chemical from those at kinetic freeze-out, in the following the former will be denoted with 
a subscript "ch.", and the latter with the subscript "f.o.", as above. 

Given a set of hadronic masses rrii and decay widths rV+j^..., the extracted values of 
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Figure 1. Values for temperature T and chemical potential /!# extracted from particle 
ratios. Different collision energies and systems are shown in different colors. Different 
symbols indicate differences in the fitting procedure, for more details see listed references 
(numbering not identical with that in main text). 



(T, Hb)c\i. displayed in Fig. |TJ are relatively insensitive to errors in the experimental deter- 
mination of total particle numbers. This can be seen in the non-relativistic approximation, 

ni (T,^)~exp(^-^) , (20) 



T 

such that Eq. ( |T8"D becomes 
Ni ( yu^ch. - yUj,ch. rrii - rrij 



n- ~ exp I— 5£ T^r) ■ (2l) 

Errors in the determination of Ni and Nj influence the fitted values for T c h. and /^s, c h. 
only logarithmically. 

On the other hand, as can be also seen from Eq. (0), the values for the hadronic masses 
rrti entering the fit influence T c h. and /j,b, ch. linearly. Most fits assume that the hadronic 
masses and decay widths at chemical freeze-out are the same as in vacuum. At such 
large temperatures and chemical potentials, however, a change of mass and decay width 
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due to in-medium interactions cannot be excluded p4|-p6|1 , m 8 — > m*, T^jk — > 
These interactions are elastic collisions which still occur at chemical freeze-out, since the 



system remains kinetically equilibrated. In ||25|| , such an analysis was performed for PbPb 
collisions at CERN-SPS energies in the framework of a chiral model. In this model, 
masses in general decrease at large density and temperature, such that the value for the 
chemical freeze-out temperature is smaller than for a fit with vacuum masses, T c h. — 144 
MeV. Interestingly, in the calculation of [25[] it turns out to be approximately the same as 
the kinetic freeze-out temperature, although at kinetic freeze-out, particles have to attain 
their vacuum masses, due to the absence of any kind of interaction. Clearly, more work 
has to be done to understand and possibly refine the freeze-out picture, for instance, also 
to include in-medium modifications of the decay widths. 



5. CONCLUSIONS 

In this talk, I first discussed single-inclusive, invariant particle spectra in a 2 — > N 
process. I explicitly demonstrated that the momentum space volume of the final-state 
particles gives rise to the exponential nature of these spectra. For the special case of 
massless particles, I showed that the slope of the spectra differs from the one for a ther- 
mally equilibrated system of particles, i.e., that exponential particle spectra alone are not 
indicative for the existence of thermal equilibrium. 

I then argued that the slope characteristics for different particle species suggest col- 
lective motion in AA collisions. Single-inclusive, invariant particle spectra then yield 
information about the average temperature Tf. Q . and average collective flow velocity (vt) 
at kinetic freeze-out. 

Finally, I discussed the extraction of the chemical freeze-out parameters T c h. and fiB,ch. 
from particle ratios. This analysis is based upon rather restrictive assumptions, and 
depends sensitively on the hadronic masses and decay widths. 

Nevertheless, it is astonishing that a fit with essentially two parameters can reproduce 
a multitude of data with reasonably good quality. There could be several explanations. 
First, the state preceding chemical freeze-out is indeed in thermodynamic equilibrium. 
The degrees of freedom could be hadrons, but since T ch is close to the temperature of 
the confinement- deconfinement transition, it is, however, also conceivable that at some 
earlier stage in the evolution of the system a thermodynamically equilibrated QGP has 
existed |TRE7|. From thermodynamic arguments alone, however, one can never distinguish 



between these two possibilities: by definition, a state in thermodynamic equilibrium has 
no knowledge about its past. 

The discussion of section ^] allows for a second explanation. Multiparticle production 
saturates the available phase space, such that final-state hadrons have exponential spectra 
with a slope that can be interpreted as a temperature, and an absolute normalization that 
can be interpreted as a fugacity. Although the discussion of section || has shown that the 
final state need not be thermodynamically equilibrated to have these properties, it then 
at least looks like hadrons are "born into thermodynamic equilibrium" p8| . In principle, 



the chemical composition could then immediately freeze out. In this case, no additional 
assumptions about an equilibrated state preceding chemical freeze-out are necessary. As 
seen in section ||, there certainly are elastic interactions which cause collective motion, 
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but this could in principle happen after chemical freeze-out. 

To distinguish between the first and second scenario, a more thorough understanding 
of multiparticle production in AA collisions is necessary. The only way to achieve this 
is to compare AA with pp and pA collisions. It is therefore mandatory to gather more 
experimental data on the latter at comparable collision energies. 
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